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Abstract 

A numerical scheme is developed for solution of the Goursat prob- 
lem for a class of nonlinear hyperbolic systems with an arbitrary 
number of independent variables. Convergence results are proved for 
this difference scheme. These results are applied to hyperbolic sys- 
tems of differential-geometric origin, like the sine-Gordon equation 
describing the surfaces of the constant negative Gaussian curvature 
(-fC-surfaces) . In particular, we prove the convergence of discrete 
-fC-surfaces and their Backlund transformations to their continuous 
counterparts. This puts on a firm basis the generally accepted belief 
(which however remained unproved untill this work) that the classical 
differential geometry of integrable classes of surfaces and the classi- 
cal theory of transformations of such surfaces may be obtained from 
a unifying multi-dimensional discrete theory by a refinement of the 
coordinate mesh-size in some of the directions. 
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1 Introduction 



The development of the classical differential geometry led to introduction 
and studying various classes of surfaces which are of interest both for the 
internal differential-geometric reasons and for application in other sciences. 
We mention here minimal surfaces, constant curvature surfaces, isothermic 
surfaces, to name just a few general classes. A rich theory of such surface 
classes is, to a large extent, a classical heritage; more recently, also numerical 
methods appeared, stimulated by applications in sciences and in scientific 
computation, including visualization. Mostly, these numerical methods are 
of a variational nature, applicable to problems described by elliptic partial 
differential equations, like the Plateau problem in the theory of minimal 
surfaces (see, e.g., PPI, ffil, pH|). It seems that no general numerical methods 



have been developed for objects of the differential geometry described by 
hyperbolic differential equations, like the surfaces with constant negative 
Gaussian curvature. 




Figure 1: Surfaces and their transformations as a limit of multidimensional 
lattices 



The characteristic property of various special classes of surfaces studied 
by the classical differential geometry turns out to be their integr ability, in 
the sense of the modern theory of solitons. One of the manifestations of 
integrability is the existence of a rich transformations theory, sometimes uni- 
fied under the names of the Darboux-Backlund transformations. Classically, 
the theory of surfaces and the theory of their transformations were dealt 
with separately to a large extent. Recently, it became clear that both the- 
ories can be unified in the framework of the discrete differential geometry 
(cf. |jS|, |BP2|| ). In this framework, multidimensional lattices with certain 



geometrical properties become the basic mathematical structures. As the 
lattice becomes more and more dense in some of the coordinates directions 
(the mesh size e — > 0), it approximates the smooth surface. The directions, 
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where the mesh size remains constant, correspond to the transformations of 
smooth surfaces (see Fig. jT]). 

The discrete differential geometry is nowadays a flourishing area which 
parallels to a large extent its classical (continuous) counterpart. Many im- 
portant classes of surfaces have been discretized up to now, see a review 
[BP2|| . Their properties are well understood. Fig.^ [| shows an example 



m 



of a continuous Amsler surface (which is a surface with constant negative 
Gaussian curvature) and its discrete analog. The characteristic property of 
discrete surfaces F : (eZ)^ i— >• with constant negative Gaussian curvature 
is that for each (x, ?/) e (eZ)^ the five points F{x,y) and F{x ± e,y ± e) lie 
in a plane. The subclass of Amsler surfaces is singled out by the condition 
that the surface should contain two straight lines. 




Figure 2: A countinous and a discrete Amsler surfaces 



Considering pictures like this, one is faced with a striking qualitative sim- 
ilarity of the continuous surfaces and their discrete counterparts. Moreover, 
from numerical experiments it became clear that the approximation is also 
quantitative: the points on the discrete surface converge to the corresponding 

"'^We are thankful to Tim HofFmann for producing these figures. 
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points on the continous surface as the mesh size e goes to zero. The typical 
picture of the approximation error is shown on Fig. |^. 




Figure 3: Maximum norm of the error F'^ — F vs. mesh size e. F : [0,r] x 
[0,r] is a smooth K-surface, : [0,r]^ x [0,r]^ ^ are discrete K- 

surfaces. (See section ^ for notations.) The lower line corresponds to r = 1.0, 
the upper one to r = 4.0. 



For a "sufficiently typical" family of discrete surfaces, we have plotted 
their maximal point distances from the respective continous counterpart. 
(See section || for details.) The slope of both curves is very close to one, 
which indicates linear convergence. 

All this suggests that it might be possible to develop the classical differ- 
ential geometry, including both the theory of surfaces and of their transfor- 
mations, as a limit of the discrete constructions, just by refining the mesh 
size in some directions. This is a common belief now, having however the 
status of folklore only, since there are no rigorous mathematical statements 
supporting it. On the other hand, the good quantitative properties of ap- 
proximations delivered by the discrete differential geometry suggest that they 
might be put at the basis of the practical numerical algorithms for compu- 
tations in the differential geometry. Again, absence of mathematical results 
on the quality of approximation prevents one from doing this. 

The present paper aims at closing this gap. We take a step in this di- 
rection, developing a numerical scheme for a class of nonlinear hyperbolic 
equations, and proving general results on its convergence. It should be said 
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that we develop the numerics, having in mind also the above mentioned 
theoretical applications to an alternative foundation of the differential ge- 
ometry. That means, first, that the class of hyperbolic systems we consider 
here includes those coming from geometric applications, and, second, that 
our discretizations respect the geometric structures, i.e. belong to the field 
of discrete differential geometry. In particular, the notion of integrability, in 
the guise of the multi-dimensional compatibility of discrete equations ||BS|| , 
will play a considerable role in our approach. 

Example. We shall illustrate our constructions by the well-known Sine- 
Gordon equation: 

5x^3/0 = sin 0. (1) 

A naive discretization of the Sine-Gordon equation could be obtained from 
([I|) by replacing partial derivatives by their difference analogs: 

5^5^0 = sin0, (2) 

where we introduced the following notation: 

^Ip{x^ y) = - {p{x + e,y)- p{x, y)^ , y) = - (^p{x, y + e) - p{x, y)^ . 

(3) 

In length, Eq. (|]) reads: 

0(a; + e,y + e) - (l){x + e, y) - (f){x, y + e) + (f){x, y) = sin 0(a;, y). 

We shall prove an approximation theorem which implies that, on finite do- 
mains, the solutions of a Goursat problem for (^) converge with e — > to 
the solutions of a Goursat problem for ([1|), provided the initial data on the 
characteristic lines converge. However, this is not the whole story. Solutions 
of the Sine-Gordon equation correspond to surfaces with constant negative 
Gaussian curvature. The discretization is non-geometric, and it not clear 
how to construct discrete surfaces from its solutions. This is closely related 
to the fact that this discretization does not inherit the integrability of the 
Sine-Gordon equation. There exists a different one, due to Hirota Q, which 
is itself a discrete integrable system. Its geometric meaning and relation to 
discrete surfaces with the constant negative Gaussian curvature was clarified 
by |BP1| ]. The Hirota's discretization of the Sine-Gordon equation reads: 



sin + e,y + e) - (j){x + e, y) - y + e) + y)) 

1 

= -^sm-{4){x + e,y + e) + (j){x + e,y) + 4){x,y + e) + 4){x,y)). (4) 
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Our theory yields the convergence results for the solutions of this equation 
to the solutions of (|lD, which can be extended to the convergence of surfaces, 
their associated families, and their Backlund transformations. The existence 
of Backlund transformations is considered as a characteristic property of 
integrability ||RS| . 

It should be noticed that there exists an extensive literature dealing with 
the numerical solution of the sine-Gordon and similar equations, see, e.g. 
[AHS1| , |AHS2| , [l^'VI , [F5| , |SV| ]. Most of these references have a physical back- 
ground and motivation, and due to this the problems settled and solved there 
are different from those relevant to the differential geometry. In particular, 
neither reference contains convergence results for the Goursat problem. No- 
tice, further, that our results are applicable to a large class of hyperbolic 
systems describing various further geometries. 

The structure of the paper is the following. In Sect. ^ we formulate the 
continuous and discrete setup of the two-dimensional hyperbolic systems and 
the corresponding Goursat problems. The C^-convergence result is proven 
which holds for all difference schemes with a local approximation property. 
The C^-approximation under the appropriate conditions is established in 
Sect. 0. The theory is extended to the case of three independent variables in 
Sect. ^. At this point the notion of three-dimensional compatibility starts 
to play the key role; it turns out to be intimately related to the integrability. 
Therefore, the convergence result holds only for difference schemes with these 
properties. The theory is illustrated in Sect. ^, where we apply the conver- 
gence results to an integrable discretization of the sine-Gordon equation, and 
thus prove the convergence of discrete i^-surfaces and their Backlund trans- 
formations to the continuous counterparts. Finally, in Sect. ^ the theory is 
extended to the case of an arbitrary number of independent variables. The 
Appendix (Sect. |^) contains the technical proofs of some statements of the 
main text. 

Further pictures of discrete K-surfaces (like the ones that constitute the 
data for Fig.^), as well as a movie visualization of the convergence can be 
found on 



http : / /www-sf b288 . math . tu-berl in . de/ ^bobenko 



2 Two-dimensional theory 

In this section we prove an approximation theorem for a certain class of 
hyperbolic differential and difference equations in two dimensions. Later on 
we will consider also more general rf-dimensional systems. The notations 
we will use for domains of independent variables are the following: let r = 
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(ri, . . . , Trf) consist of positive numbers r-i > 0, then 

n(r) = [0, n] X ... X [0, u] C R'. (5) 

As domains for discrete equations, we use parts of rectangular lattices inside 
f^(r), with possibly different grid sizes along different coordinate axes e = 
(ei, . . . , erf): 

d 

Q^{r) = [0, X ... X [0, r.Y'^ C l[{e,Z). (6) 

i=l 

where [0,rY = [0, r] fl (eZ). The dependent variables of the differential and 
difference equations under consideration are supposed to belong to a normed 
vector space X, the norm in which will be denoted by | ■ |. 

In the two-dimensional situation the notations will be, somewhat incon- 
sistently, simplified: we denote continuous domains by Q{r) = [0,r] x [0,r] C 
R^, and discrete ones by f2^(r) = [0,rY x [0,r]^ C (eZ)^, Each fi^(r) contains 
0(e^^) of grid points. If 62 is an integer multiple of ei, then Q'^^{r) C Q'^^{r). 
It will be convenient to assume that e attain only values of the form 2"^^ 
with a positive integer k; for such e, the relation ei < €2 implies that €2 is an 
integer multiple of ei. We also define the limiting domains 

fi°(r) = U fi^(r), (7) 

which are everywhere dense in Q{r). Each point x G Q^{r) belongs to Q^{r) 
with e = 2~'^ for all k large enough, A; > ko{x). This allows us to speak about 
convergence by e = ^ of sequences of functions a*^ defined on Q'^{r). In 
case of convergence, the limiting function may be thought of as defined on 
Q^{r). If such a limiting function is continuous (Lipschitz) on Q^{r), it can 
be extended to a continuous (resp. Lipschitz) function on Q{r). 

Introduce the difference quotient operators 6^ and 6y, acting on functions 
p : fi^(r) X, by the formulas (||). The functions 6lp and 6yP are defined 
everywhere on fi^(r) except for the points with maximal x or y coordinate. 

Definition 1 By a continuous 2D hyperbolic system we call a system 
of partial differential equations for functions a, b : Q{r) ^ X of the form 

d^a = f{a,b), dyb = g{a,b), (8) 

with smooth enough functions f,g:XxX—>-X. A Goursat problem 

consists of prescribing the initial values 

a(x,0) = ao(x), b{0,y) = boiy) (9) 

forx G [0,r] and y G [0,r], respectively. The functions aQ,bQ : [0,r] i-^ X are 
also supposed to be smooth enough. 
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Definition 2 A discrete 2D hyperbolic system (or, better, a one- 
parameter family of such systems ) consists of two partial difference equations 
for a, b : Q'^{r) X of the form 

6la = f^{a,b), 6lb = g^{a,b), (10) 

with smooth functions f^, g"^ : X x A! X . A Goursat problem for this 
system consists of prescribing the initial values 

a(x,0) = a^(x), b{Q,y) = bl{y) (11) 

for X G [Ojr]'^ and y G [0,r]'^, respectively. 

Example. Any equation of tlie type dxdyU = F{u) can be brought into 
tlie form (H) by a variety of substitutions a = A{u,dxu) and b = B{u,dyu). 
For instance, a canonical way to do this for the Sine-Gordon equation ([l|) is 
to introduce two new dependent variables 

a = dx(j), b = (j), (12) 

which have to satisfy a system of the form (|^): 

dya = smb, d^b = a. (13) 

The second order difference equation (^ can be dealt with in a way which 
mimics the continuous situation: introduce two new dependent variables 

a = 61<P, b = 0, (14) 

then they have to satisfy a discrete 2D hyperbolic system 

Sya = smb, 5lb = a. (15) 

A similar procedure can be performed with the discretization (H). Set 



5^0, b = <p+'-6l<p. (16) 



(The second equation may be written as 6(x, y) = (0(x, y + e) + (j){x, y))/2.) 
Then the following holds: 

b{x + e,y)-b{x,y) = ^{a{x,y + e) + a{x,y)), (17) 



iea{x,y + e)/2 _ ^iea{x,y)/2 ^ ^ ^^ib{x + e,y) _ ^-ib{x,y)y (^^g^ 



(the first equation guarantees the existence of (f) for a given pair of functions a, 
6, satisfying (|T6|), while the second one is equivalent to @). These equations 
can be solved for a(x, y + e) and h{x + e,y). The result reads: 

2, 1 - (eV4) exp(-i6 - iea/2) e„ 

Both discrete 2D hyperbolic systems (|15|), (|19|) approximate the continuous 
one fllSD in the sense of the next definition. 



Definition 3 A discrete 2D hyperbolic system approximates the con- 
tinuous one if the functions f^, g"^ continuously depend on e E (0,eo), 
and 

fia, b) = hm f (a, 6), gia, b) = \iu,g%a, b), (20) 

uniformly on any compact subset of X x X. If this convergence holds in the 
-topology, we speak about the approximation. If 

r (a, b) = fia, b) + 0(e), g^ia, b) = g{a, b) + (9(e) (21) 

uniformly on any compact subset of X x X , we speak about the 0{e)- 
approximation. 

The following result is almost obvious: 

Proposition 1 The Goursat problem for a discrete 2D hyperbolic system 
( fZT^ has a unique solution [a'^^W) on fi^(r). 

Proof. At this point a useful remark should be done. Although our 
notations might suggest that the variables (a, b) are attached to the points 
of the two-dimensional lattice f^'^(r), it is more natural to assume that they 
are attached to the edges of this lattice: a{x,y) - to the horizontal edge 
connecting the vertices {x,y) and {x + e,y), and b{x,y) - to the vertical 
edge connecting the vertices {x, y) and (x, y + e). So, the equations ([T0| ) give 
the fields sitting on the right and on the top edges of an elementary square, 
provided the fields sitting on the left and on the bottom ones are known. 
See Fig. ^. By induction, the whole solution can be calculated, starting 
with the fields sitting on all edges on the coordinate axes. ■ 

Now the main result of this section can be formulated. 

Theorem 1 Let a family of discrete 2D hyperbolic systems ( |7Z[ ) 0{t)- 
approximate the continuous 2D hyperbolic system ^ in the sense. Let 
also the discrete initial data ( [71[ j approximate the continuous ones 

al{x) = ao{x) + 0(e), bl{y) = bo{y) + 0{e) (22) 
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(x, y + e) (x + e, ?/ + e) 



a{x,y + e) 



h{x,y) 



b{x + e,y) 



a{x,y) 




+ e,y) 



Figure 4: An elementary quadrilateral 



uniformly for x G [0,r]^ and y G [0,r]^, respectively. Then the sequence 
of solutions (a*^, 6^) converges uniformly in the following sense: there exist 
f G (0,r] and a pair of Lips chitz- continuous functions {a,b) on Q{f) such 
that 



for all {x, y) G fi^(f) . The functions (a, b) solve the continuous Goursat prob- 
lem for ^ on VL{f). If the functions f , g admit global Lipschitz constants, 
then one can choose f = r. 

Proof. In general one cannot expect f = r in the Theorem because the 
solutions of the limiting equations may have blow-ups that are absent in the 
discretization. Consequently, one essential step in proving the theorem is to 
attain e-independent a priori bounds on a*^ and 6"^. 

Lemma 1 Let the norms of initial data a^, 6g be bounded by e-independent 
constants. Then there exists f G (0, r] such that the norms of the solutions 
(0^,6*^) are bounded on the respective Q'^{f ) independently of e. Furthermore, 
if f and g have a global Lipschitz constant C on the whole of X x X , then 
one can choose f = r. 

Proof of Lemma |I|. Let IoqI, |&ol ^ ^0 with Mq > 0. We show that, fixing 
an arbitrary Pq > Mq, it is possible to find f < r (independent of e) such that 
one has |a^(x,?/)|, |6^(x,?/)| < Pq for {x,y) G Q''{f). Actually, it is enough to 





(23) 



take 



r = 



Pq-Mq 



(24) 



where 



.F(M) = sup sup {\f^{a,b)\,\g'{a,b)\}. 



(25) 



e \a\,\b\<M 
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Indeed, from the difference equations (|T0), written in length as 



a'{x,y) = a'{x,y-e) + ef'{a'{x,y-e),b'{x,y-e)), (26) 
b'{x,y) = ¥{x-e,y) + eg'{a'{x-e,y),¥{x-e,y)), (27) 

we can conclude by induction that 

\a'ix,y)\<Mo + yJ'iPo), |6^(x, y)| < Mq + x^(Po), (28) 

at least as long as the left-hand sides remain < Pq. Since Mo + fjF(Po) = Po, 
they remain < Pq for all x,y < f. 

The argument for the case when /, g, and therefore f^, g'^ have a global 
Lipschitz constant C, is a little bit different. Set 

A{x,y) = max{\a'{x,y)l\¥{x,y)lMo}. (29) 

From (p6[), (pTf ) it follows readily that 

A(x, y) < (l + eJ(A(x', y'))) ■ A{x' , y'), (30) 

where {x',y') = {x,y — e) or {x',y') = (x — e, ?/), depending on which of 
\a''{x,y) \ or \b''{x,y)\ is greater, and where 

J(M) = i- sup sup {|f(a,6)|,|/(a,6)|}. (31) 

Now, this function admits for M > Mq an estimate by an absolute constant: 
J{M) < ^(ir(0,0)| + 1/(0,0)1) + 2/: := /C. 

Hence, 

A(x,t/)<(l + e/C)-A(x',i/')- 
Using Lemma below (or a simple induction), we find: 

A{x,y) < Moexp(2/C(x + y)) < Moexp(4;Cr). (32) 

This finishes the proof. ■ 

On the last step we used the simple particular case = 2, k = of the 
following lemma, which will be used repeatedly later on. Its proof is given 
in the Appendix. 
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Lemma 2 Let the function A : f2^(r) M+ satisfy the following condition: 
for any x G r2^(r) different from the origin, there exists an index 1 < i < d 
such that 

A{x) < {1 + eilC)A{x - eiei) + eiK, (33) 

where JC, k are some non-negative numbers, and is the unit vector of the 
i-th coordinate axis. Then 

d 

A(x) < max (^A(O), exp (^2/C^xj). (34) 

We demonstrate further that the solutions (a'^, 6^) are not only uniformly 
bounded, but are actually Lipschitz continuous with a Lipschitz constant 
independent of e. In virtue of the equations ([l0[) and Lemma [l| it is clear 
that the difference quotients dytt" and S^ff are uniformly bounded. This turns 
out to be true also for 5ia'' and 55, 

Lemma 3 Let the initial data of the continuous Goursat problem ao? • 
[0,r] ^ X be functions, with the C^-norm less then M > 0, and let the 
initial data of the discrete Goursat problem a^^h^ : [0,r]^ i— > X satisfy 

\al{x) - aoix)\ < Me, \bl{y) - boiy)\ < Me. (35) 

With f G (0, r] chosen according to Lemma the expressions d^.a'^ and 5yb^ 
are bounded on the fi'^(f) by e-independent constants. 

Proof. Let Mi be a common bound on the values of the solutions of the 
discrete Goursat problems. Set 

M2 = sup sup {\na,b)\,ma,b)\,\dana,b)\,...,\dtg'{a,b)\] (36) 

e |a|,|6|<Afi ^ 

(recall that f^ ^ f and g'^ ^ g locally uniformly in C^). One can assume 
that M from the condition of the Lemma is greater than Mi and M2. By 
the mean value theorem, we find an estimate for 6l.a'^{x,y) with y = 0: 

|5^ao(x)| < \6lao{x) \ + e~"^|ao(x + e) — ao{x + e)| + e""^|ao(x) — al{x)\ < 3M. 

(37) 

Proceeding from y to y + e, we find: 

\6la'{x,y + e)\ < \6la'{x,y)\ + e\6ir{a'{x,y),¥{x,y))\ 

< \6la%x,y)\ + eM{\Sla%x,y)\ + \6l¥{x,y)\) 

< {l + eM)\6la'{x,y)\+eM^. 
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Now Lemma |^ yields the desired estimate: 

\Sl^a'{x,y)\ < 4Mexp(Mf). 
The same reasoning apphes to 5^6". ■ 

Proof of Theorem |T], continued. We have to show that the e- 
dependent solutions of difference equations have limits (a", W) on Q^{f) with 
the convergence rates 

sup \a'ix,y) - a'^ix,y)\,\¥ix,y) -b%x,y)\ = 0{e), (38) 

and that these limiting functions (a°, 6") can be extended to continuous func- 
tions (a, b) on Q{f) satisfying the differential equations. 

Take M > M2 (cf. (0)) such that it bounds also a"^ and 6'^ along with their 
respective difference quotients (see Lemma H). We will prove that a''{x,y) 
and h^i^x^y) are Cauchy sequences at any point (x, y) G Vf{f). To this end, 
fix e, e' which are of the usual form 2"^^; we assume e > e', so that e/e' is a 
natural number. Below we use the following abbreviations: A^^ a (a;, y) for 
|a^(x, ?/) —a'' (x, ?/)|, and /"^[a"^, fe'^](a;, ?/) for f''{a'^{x, y), lf[x, y)), etc. We have: 

A"'a(x,i/) < 

A«'a(x,y - e) + erK,6^](x,y - e) - e' ^ r'[a^ &^'](x, y - «:e') 
A"'a(x, y - e) + e' ^ I [aS 6^] (x, y - e) - [a^' ,¥']{x,y - ^e') | . (39) 

K = l 

For the terms in the last sum with a fixed < k < e/e' we have: 

\r[a%¥]ix,y- e) - f [a'\¥']ix,y - Ke')\ < 

\r[a%¥]ix,y- e) - r'[a^,¥]ix,y - e)| (40) 

+ \f'[a^,b^]{x,y-e)-r'[a^',b^']{x,y-e)\ (41) 

+ |r'[a^',6^'](x,y - e) - r'[a^',¥']{x,y- Ke')\. (42) 

Here = 0{e) because f = f + 0{e), (||) = C(e) due to Lemma |, and 
dg) < M (^A"'a(x, y-e) + A"'6(x, y - e; 



Putting all this on the right-hand side of (|39|), and taking into account that 
the sum there contains 0{e/e') terms, we find: 

A"'a(x, y)<{l + eM)A"'a(x, y - e) + eMA"'6(x, y - e) + 0{e^). 
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The analogous estimate holds for b: 

A"'6(x, < (1 + eM)A"'6(x - e,y) + eMA"'a(x - e,y) + 0{e^). 
Introducing the quantities 

mo := sup {|ao(-) -ao(-)U&o(-) -&o(-)l K 

A(x, := max{A"'a(x, y), A"'6(x, y), mo}, (43) 

we get the estimate which holds also for the points of the boundary of Q'^{f) 
different from the origin: 

A(x, y)<il+ eM)A{x', y') + ^(e^), (44) 

where either {x',y') = {x,y — e), or {x',y') = {x — e,y). Now Lemma ^ 
yields the final estimate: A{x,y) = 0{e). Thus a''{x,y) and b''{x,y) are 
proved to form Cauchy sequences at every point {x,y) G Q^{f). Their limits 
a°(x, y) and 6°(x, y) satisfy (|38[). These limiting functions obviously have the 
same Lipschitz constants as all the a*^ and 6^. Therefore there is a unique 
(Lipschitz) continuous extension of these functions from f2°(f) to f2(f); we 
call these extensions a and b. 

It remains to be shown that (a, b) solve the differential equations (|^). To 
do this, we prove that certain integral equations hold. Let {x,y) G Q^{f), 
then: 

s/A-i 

a'{x,y) = al{x) +e ^ /'[a', 6'](x, /ce) = 

A:=0 

y/e-1 y/e-1 

al{x)+e J2 f[a,b]{x,ke)+eJ2 (r [aS 6^](x, A;e) - /[a, 6](x, A;e)) . 

k=0 k=0 

All the terms inside the last sum are uniformly 0{e), therefore the whole 
sum (with the pre-factor e) is also 0{e). The first sum on the right-hand 
side is, up to 0{e), equal to f{a{x,ri),b{x,ri))dri, which exists since / is 
smooth and a, b are Lipschitz. Therefore, 

ry 

a{x,y) = ao{x) + f[a,b]{x,r])dr] (45) 
Jo 

on f2°(f). By continuity, this holds on all of Q{f). It follows that a is 
everywhere different iable with respect to x, and d^a = f{a,b). The function 
b is treated in the same manner. ■ 
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Corollary 1 The two dimensional hyperbolic Goursat problem pos- 
sesses a unique classical solution. 

Proof. The existence part of this statement is an immediate consequence 
of Theorem ^ uniqueness is easy to show as follows. Let (a, h) and (a, h) be 
two solutions of the system of integral equations consisting of (|^) and 

h{x,y) = h,{y)+ [ g[a,b]{^,y)dt (46) 

Subtracting the corresponding equations, we find after some simple manipu- 
lations: 

A{x,y)<LQ^ A{^,y)d^ + J\{x,r])dr]y 

where we introduced the deviation function 

A{x,y) = \a{x,y) - d{x,y)\ + \b{x,y) - b{x,y)\, 

and L is a common Lipschitz constant of /, g over the range of values of 
a, b. Now the 2D version of the classical Gronwall inequality {d = 2 case of 
Lemma ^ below) implies A{x,y) = 0. ■ 

Lemma 4 Let the continuous function A : Q{r) satisfy 

A(xi, . . . ,Xd) < / A{xi, Xj_i, ^j, Xj+i, Xd)dij + Q, (47) 

with some constants L,Q >0. Then 

d 

A{xi,...,Xd) <2Qe^^{2dLY^x^. (48) 



1=1 



Proof of this lemma is put in the Appendix. 

3 Additional Smoothness 

In this section we show that the discrete techniques can be used to prove 
regularity of the solutions to the hyperbohc equations. 
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Theorem 2 Let the assumptions of Theorem hold; additionally, let the 
partial difference quotients up to the order k + 1 of the discrete initial data 
be uniformly bounded independently of e: 

mral{x)\, \{Sirbl{y)\ < M, m,n<k + l, (49) 

Suppose that the convergence f^—^f,g^-^g is locally uniform in C^^^ . 
Then the limit functions a = lime_>o a^, b = liiiie^o^^ belong to , and 
moreover these limits are uniform in : 

sup m^{5lTa^-d^d^a\A{5ir{5irh'-d^d^b\-.Q, m + n<k. 

(50) 

Remark. Assume that the convergence f'^—>-f,g'^—>-g is locally 
uniform in C°° and that the initial data ao,&o are C°°-smooth. Then the 
canonical choice of the discrete initial data aQ^x) = ao{x) and &o(2/) — ^oiv) 
guarantees the C'^-convergence for any k. One may then loosely speak of 
C°°-approximation. 

The ideas of the proof and even the essential estimates are basically the 
same as one would most likely use in the continuous setting. However, the 
discrete setting has the advantage that in contrast to higher-order partial 
derivatives all the difference quotients automatically exist, commute with 
each other etc. 

First, one obtains a priori estimates not only for the values of a"^, 6*^, but 
for their higher order difference quotients. We will need discrete analogues of 
the C'^-norms. Let 3^ be a normed linear space; for a function u : fi'^(r) y, 
define 

||m||o = sup (51) 

and 

The following statement comes to replace the chain rule, which is no more 
available in the discrete context. 

Lemma 5 Let f : X ^ X ^ X be a smooth function, and consider two 
functions a,b : Q'^{r) — > X. Then the K-th order difference quotients of 
f{a, b) can be estimated as follows (m + n = K): 

m)"'i5irf[a,b]i^,y)\ < A- (|(5^)™(5^)"a(x,i/)| + mri5irbix,y)\) + 

Pi\\a\\K-iA\b\\K-i) + Be. (53) 
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Here P is a polynomial on two variables of total degree < K with positive 
coefficients; the constants A, B, and the coefficients of P depend only on K, 
on the Lipschitz constant of the functions a and h, and on the expressions 
\\{D^M)Mhfork + l<K + l. 

The technical proof of the lemma in put in the Appendix. 

Remark. In applications of Lemma ^, we have to handle a whole e- 
dependent family of functions a'',b'', g" at once. Thus the constants A, 
B and the polynomial P become, in principle, e-dependent, too. However, 
under the conditions like those of Theorem H, one can choose these constants 
independently of e. In particular, if sup^ ll^^'^lli^' ^-^d sup^ ll^'^llx are finite, then 
{5l.)'^{5yY f'i^a'^ , ¥) is bounded for m + n = i^', independently of e. 

Lemma 6 Under the conditions of Theorem 

sup ||a''||fc_|_i < oo, sup < oo. (54) 

Proof goes by induction with respect to the total degree m + n = K<k + l 
of the difference quotient. So assume that sup^ ||a''||i4'_i, sup^ < oo is 

already proved. 

First, look at the difference quotients of a*^ with n > 0. Then 

{5ir{5ira\x,y) = {5ir{5lT-'f\a\¥){x,y), (55) 

and the statement follows from the induction hypothesis by using Lemma |^, 
as pointed out in the preceeding remark. Analogously, the difference quo- 
tients of h'^ with m > are readily estimated. If n = 0, then we have: 

{5ira^{x, y) = (5^)"^a^(x, y - e) + e(5^)"^r [a^ ¥]{x, y-e). (56) 

Apply Lemma ^ to find: 

\{Sira%x,y)\ < il + eAmra'{x,y-e)\ 
+eA\{6:,r¥{x,y)\ 
+eP{\\a'\\K-i,m\K-i) + e^B. 

By the induction hypotheses, the case considered above and the remark 
following Lemma ^, the quantities on the right-hand side are controlled, 
so that Lemma ^ can be applied, which yields a uniform, e-independent 
estimate for ((5^)'"a^ Completely the same reasoning applies to (5y)"6^ ■ 

Theorem ^ is now a direct consequence of the next lemma. 
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Lemma 7 Let a sequence of functions {u" : fl''{r) X}^(zEq be bounded in 
the discrete C^^^-norm, independently of e: 

sup llu^llfc+i = M < oo. (57) 

Then there exists a function u G C''{Q{r),X), and a subsequence {u^jeg^;^ 
for which 

snp\{6ir{6iru'{x,y)-d:^d;;u{x,y)\^0 as e 0, e e E^. (58) 

Proof. The set Q^{r) consists of countably many points; choose any enu- 
meration 

l]°(r) = {(x„,y„)|n = 1,2,3,...}. (59) 

Since all difference quotients up to the order k are bounded at {xi,yi), there 
exists a subsequence {w'^jee^Bi such that 

i5lriSlru^{x^,y^)^u^^'^\x^,y^) as e ^ 0, e e E^. (60) 

for all m + n < k. Repeat this procedure at {x2,y2), where a convergent 
subsubsequence with E2 C Ei is selected, and so on. We get a series of 
infinite sets -Ej+i C E^. Eventually, for all m + n < the limits u^"^'^\x, y) 
are defined everywhere on Vf[r). Since < M holds e- uniformly, the 

functions thus defined have the Lipschitz property on f2*^(r), so they 

possess unique continuous extensions to n(r) with M as a Lipschitz constant. 

We are ready to describe the set E^o from assertion of the lemma. We 
construct it as an infinite sequence of numbers converging to zero, in such 
a way that 

sup \m^{5l^Tu\x,y)-u'^^^^\x,y)\ < 2-^(2M + 1). (61) 

Let eo = 1, and take < ej_i of the usual form (integer power of 1/2), 
satisfying 

\{5lT{5yru^^{x\y')-u^^^''^\x\y')\ < (62) 

for all m + n < A; at all points (a;', y') G Vt^ ' (r) (we have with necessity 
ej < 1~\ therefore ^ ^ ir) C f2'^j (r)). Such exists, because for finitely 
many lattice sites of ' there holds a pointwise convergence argument 
from the beginning of this proof. In order to establish the estimate (^) on 
the larger set n'^^(r), we make use of the Lipschitz properties. Let [x^y) G 
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fi'^-'(r) be arbitrary; then there always exists a point {x',y') G fi^ ^(r) with 
\x' — x\ + \y' — y\ < 2~K Thus we conclude 

+ \{m^{5l^ru^^{x\y')~u^-'-\x\y')\ 
+ \u'^'^^''\x\y')-u^'^^''\x,y)\ 
< M2-~^ + 2-^ + M2-K 

It remains to show that are indeed the respective partial derivatives of 

u. This can be done by proving the corresponding integral equations: from 

x/t-l 

(SinSiru^x, y) = {STi^iru^O, y) + ^J2 {^ir^^^lTu^ike, y) (63) 

fc=0 

there follows in the limit e — > 0, e G E^o 

= n('"'")(0,i/) + r u^'^^''-\U)di, (64) 

^0 

so that (9^m(™'")(x,?/) = m(™+1'")(x,?/). ■ 

Proof of Theorem 0. With Lemma |^ at hand, the convergence of 
the difference quotients of a^, 6'^ is proved as follows. Any subsequence of 
(a*^, has a subsubsequence which converges uniformly with all its difference 
quotients to a limit function. But by Theorem |l], this limit function is 
the unique solution (a, h) of the correspondent continuous Goursat problem. 
If any subsequence of a given sequence has a subsubsequence converging to 
one and the same limit, then the whole sequence must converge to that same 
limit. ■ 



4 Three— dimensional theory: 

approximating Backlund transformations 

The Sine-Gordon equation (|l]) is a very special one, in that it posseses 
Backlund transformations. This means that from a given solution we can 
construct new solutions by solving ordinary differential equations only. The 
famous formula for a family of elementary Backlund transformations i— 
for (|l|) reads: 

4>-<t) r, ~ r, 2 + ^ ^ 

dxtp + Ox(t) = 2a sm — - — , dycf) - dy(f) = - sm — - — . (65) 
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A direct calculation shows that this system is compatible (i.e. dy{dx4>) = 
dx{dy(f))), provided is a solution of the Sine-Gordon equation, and then 
is also a solution. An equivalent way to express this state of affairs is to 
introduce, along with the variables a, b from (|T^, also the auxilary function 
9 = {(j) — 0)/2, which satisfies the following system of ordinary differential 
equations: 

dJ = -a + asiYie, dyO = -sm{h + e). (66) 

a 

This system is compatible in the sense that dy{dx6) = dx{dy9), provided (a, b) 
solve the system (|13D equivalent to the Sine-Gordon equation, and the initial 
value 

^(0,0) = ^^0 (67) 
defined uniquely a solution of (^). Then the formulas 

a = a + 2dx9 = -a + 2asm9, b = b + 29 (68) 

deliver a new solution (a, b) of the 2D hyperbolic system ( p!3D equivalent to the 
Sine-Gordon equation. Clearly, Backlund transformations can be iterated in 
a straightforward manner. 

Definition 4 A continuous 2D hyperbolic system with a Backlund 
tranformation is a compatible system of partial differential and difference 
equations 

dya = f{a,b), dxb = g{a,b), 

dx9 = u{a,9), dy9 = v{b,9), (69) 

6,a = aa,9), 6,b = r]{b,9) 

for functions a,b,9 : Q{r, R) X , where 

n{r, R) = {{x, y, z) \ (x, y) G n{r), ^ = 0, 1, . . . , i?}. (70) 

Here f,g,u,v,^,r] : X x X ^ X are asuumed to be smooth functions. A 
Goursat problem is posed by the requirement 

a(x,0,0) =ao(x), b{0,y,0) = bo{y), 9{0,0, z) = 9o{z) (71) 

for x e [0,r], y & [0, r], and z E {0,1, . . . , R}, respectively, with given smooth 
functions ao{x), bo{y) and a sequence 9o{0), . . . ,9o{R — 1). 
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The compatibility conditions mentioned in this definition, are to assure 
the existence of solutions of the above Goursat problem. They follow from 



In length, these conditions for ( p^D to be Backlund transformations of the 
2D hyperbolic system read: 

Dau{a, 6) ■ f{a, b) + Deu{a, 6) ■ v{b, 6) = 

= Dhv{b, 9) ■ g{a, b) + Dev{b, 6) ■ u{a, 9), 
Daiia, 6) ■ f{a, b) + De^a, 6) ■ v{b, 6) = 

= f(^a + aa,9),b + r^ib,e))-fia,b), (72) 

D,r]{b, 9) ■ g{a, b) + De7]{b, 9) ■ u{a, 9) = 

= g[a + e(a, 9),b + 7]{b, 0)) - gia, b). 

The existence of Backlund transformations may be regarded as one of the 
possible definitions of the integrability of a given 2D continuous hyperbolic 
system. For a given 2D continuous hyperbolic system with Backlund trans- 
formations, not every discretization possesses the analogous property. For 
instance, the naive discretization (Q) of the Sine-Gordon equation does not 
admit Backlund transformations, while the integrable discretization (^) does. 
The difference analogs of the formulas (^) read: 

sin + e,y)- y) + 0(x + e,y) - y)) = 

ea 1 ~ ~ 
= Ysin-(0(x + e,y) + 0(x,?/)-0(a; + e,y)-0(a;,y)). (73) 

sin ^(0(x, y + e)- 0(x, y) - (f){x, y + e) + (f){x, y)) = 

= 7^sin^(0(x,?/ + e) + 0(x,?/) + 0(x,?/ + e) + 0(x,?/)). (74) 
2a 4 



Obviously, in the limit e — these equations approximate ([65|). A very 
remarkable feature is that these equations closely resemble the original dif- 
ference equation (|), if one considers the tilde as the shift in the third z- 
direction. Upon introducing the quantity 9 = (0 — 0)/2, one rewrites (|73|), 
([7^) in the form of the system of first order equations approximating (|66D, 
(H): 

S-S . -..ilog \-<;"/!'-P(-; + -"^' . (75) 

ze 1 — {ea/2) exp{z9 — iea/2) 



xen 1 1 1 - {e/2a)exp{-ib-i9j 
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and ^ 

a = a + 2519, b = b + 29 + e5l9. (77) 

This suggests the following definition. 

Definition 5 A discrete 3D hyperbolic system is a collection of com- 
patible partial difference equations of the form 

6la = f'{a,b), 6lb = g'{a,b), 

5l9 = u'{a,9), 5l9 = v\b,9), (78) 

5,a = C{a,9), 5,b = 7]\b,9), 

for functions a,b,9 : Q'^{r, R) X , where 

n^ir, R) = {{x, y, z) \ {x, y) e n\r), z = 0,1, . . . , R]. (79) 

Here the functions f^, g'^,u'^,v'^,^'^,r]'^ : X ^ X ^ X are smooth enough. A 
Goursat problem consists of prescribing the initial data 

a{x,0,0) = al{x), b{0,y,0) = bl{y), 9{0,0, z) = 9l{z) (80) 

for X G [0, r]% y G [0, r]% and z G {0, 1, ... , R}, respectively. 

Compatibility conditions are necessary for solutions of (|75D to exist. 
These conditions express the following identities that have to be fulfilled for 
the solutions: 

51{519) = 61(619), 6l{6.a) = 5,(5^a), 6l{6.b) = 6Mb). 
In length, these formulas read: 

u' (a + ef'ia, b),9 + ev'{b, 9)^ - u'{a, 9) = 

= v'(b + eg'{a, b),9 + eu'{a, 9)^ - v'{b, 9), 
C {a + ef\a, b),9 + ev\b, 9)) - C{a, 9) = 

= er (a + r(a, 9),b + ^\b, 9)) - ef\a, b), (81) 
7]' (b + eg'{a, b),9 + eu'{a, 9)^ - r]'{b, 9) = 

= (a + r(a, 9),b + r,\b, 9)) - eg\a, b). 

This has to be satisfied identically m. a,b,9 E X . 

As demonstrated in [|B|j|| , the compatibility of a discrete 3D hyperbolic 
system is closely related to its integrability in the sense of the sohton theory. 
Moreover, such a key attribute of integrability as a discrete zero curvature 
representation with a spectral parameter can be derived from the fact of com- 
patibility. 



22 



Proposition 2 The Goursat problem for a discrete 3D hyperbolic sys- 
tem ( ff^ ) satisfying the compatibility conditions ^8l\ ) has a unique solution 
{a%¥,e') on n'{r,R). 



Proof. Like in the proof of Proposition |1|, it is enough to demonstrate 
that the solution can be propagated along an elementary "cube" of the 
three-dimensional lattice, then the induction will end the proof. Again, 
it is convenient to assume that the variables a{x,y,z), b{x,y,z), 9{x,y,z) 
are attached not to the points {x,y,z) G fl''{r,R), but rather to the edges 
[{x, y, z),{x + e, y, z)], [{x, y, z), (x, y + e, z)], [(x, y, z), (x, y,z + 1)], respec- 
tively. Denote (in this proof only) shifts of the edge variables in the directions 
ofx,y,z axes by the subscripts 1,2,3, respectively. (See Fig. |^.) Then the 
values (a2, foi) are detrmined by the first equation in (^), the values (6^1, 6*2) - 
by the second one, and the values (03, 63) - by the third one. For the each one 
of the values (023, 613, 612) sitting on the edges adjacent to {x + e,y + e, z + 1) 
we get two possible values (from two equations defined on two facets sharing 
the corresponding edge). The compatibility conditions ( pi]) guarantee that 
these two values for each of the edge variables actually coincide. ■ 




{x,y,z) a {x + e,y,z) 



Figure 5: Three-dimensional consistency 



Theorem 3 Let the family of discrete 3D hyperbolic systems ^7^ ) satisfying 
the compatibility conditions approximate the continuous 2D hyperbolic 
system with a Bdcklund transformation ([^^. Let the approximation be of the 
order 0{e) in the -sense, so that the relations 

f^{a,b) = f{a,b) + 0{e), etc. (82) 
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and similar relations for the first partial derivatives hold uniformly on any 
bounded set of {a,b,6) G X^. Let also the initial data be approximated ac- 
cording to 

a',{x) = ao{x) + 0{e), b^iy) = bo{y) + 0{e), 9^,{z) = 9o{z) + Oie) 

uniformly for x G [0, rY, y G [0, r]% z G {0, 1, . . . , R}, respectively. 

Then, for some f G (0,r], the sequence of solutions (a*^, 6*^, 6'^) has a uni- 
form limit of Lip schitz- continuous functions (a, 6, 9) on f2(r, R) in the sense 
that the relations 

a\x,y,z) = a{x,y,z) + 0{e), 
b%x,y,z) = b{x,y,z) + 0{e), 
9'{x,y,z) = 9ix,y,z) + Oie) 

hold uniformly in {x, y, z) from the respective fl^i^p, R). Furthermore, (a, b, 9) 
solve the Goursat problem for the continuous 2D system with a sequence of 
Bdcklund transformations. 

Proof parallels the proof of Theorem ^ and starts with an a priori estimate 
for a% 9\ 

Lemma 8 Let the norms of initial data a^, 6q, 9q be bounded by e- 
independent constants. Then there exists f G (0, r] such that the norms of 
the solutions {a", ¥,9'') are bounded on the respective VL^if^R) independently 
of e. If the right-hand sides of Eqs. ( jT^ possess a global Lipschitz constant 
C on the whole of X x X x X , then one can choose f = r. 

Proof of Lemma g. Let |a^|, 16*^1 < Mq with Mq > 0. Set 

^(M) = sup sup {\f^{a,b)\,...,\r]'{b,9)\}. (83) 

e |al,|6|,|6)|<M 

Choose an arbitrary Pq > Mq, and define inductively Pj+i = Pj + J^{Pj) for 
j = 0,1, . . . , R — 1. Finally, set 

Pi -Mo 

r = mm — ; — — , (84) 
j=o,i,...,R 2J^{Pj) ' ^ ^ 

so that 

Mo + 2fJ^{Pj) < Pj for all j = 0,1,...,R. (85) 
We show now by induction in z that for all {x, y, z) G VL'^{r, R) there holds: 

\a%x,y,z)\, mx,y,z)\, \9%x,y, z)\ < P,. (86) 
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Indeed, like in Lemma ^ (see (|28|)), from the first pair of equations in ([78| ) 
we conclude that 

\a'{x,y,0)\ < Mo + y^(Po), \b'{x,y,0)\ < M, + x:F{P,), (87) 

and similarly from the second pair of equations in (|78D we conclude that 

|r(a;,y,0)| <Mo + (x + 2/)^(Po), (88) 



as long as the left-hand sides of these inequalities remain < Pq. Due to (pSj), 
they remain < Pq for all (x, y) G fi'^(f). Now assuming (|86D for a given z, 
the third pair of equations in (173) immediately implies: 



|a^(x, z + 1)1, |6^(a;, ;2 + 1)| < P, + ^(P,) < P,+i, (89) 

while from the second pair we derive: 

\e\x,y,z + l)\<M^ + {x + y)T{P,+i)<P,+i. (90) 

This proves (|86|) , and therefore the first statement of lemma. The second 
one is proved similarly to the analogous statement of Lemma |I]. ■ 

Further, like in the two-dimensional case, we have to demonstrate that 
not only the solutions {a", ¥,6") are bounded on f2''(f, P), but also their 
difference quotients do. 

Lemma 9 Let the initial data of the continuous Goursat problem oq, 6o • 
[0,r] \^ X he functions, with the C^-norm less then M > 0, and let the 
initial data of the discrete Goursat problem a^.h^ : [0,r]^ i-^ X satisfy 

\al{x) - ao(x)| < Me, \hl{y) - ho{y)\ < Me. (91) 

Let f G (0, r] be chosen according to Lemma |^. Then all the difference quo- 
tients S^a", Sya", S^b^, dyh", 5^.6'% and byd'' are bounded on all of the (^"(f, P) 
by e-independent constants. 



Proof of Lemma |9|. Let the estimates (|86| ) hold for all {x, y, z) G f2'^(r, P). 
Then from the equations (|78D it is immediately seen that 



\5la% \5lh% \5ie% \5ld^\<HP.)- 

Therefore, like in Lemma |^, we only have to estimate the quotients S^a*^ and 
5'yb'^. Let M be chosen to be greater than the common Lipschitz constant of 
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the functions (7", and r^" for all {a,b,9) with \a\, \b\, \9\ < Pr. Then, 
proceeding as in the proof of Lemma we have: 

0)1 < \6la'{x, y - e, 0)1 + e\6ir[a% ¥]{x, y - e,0)\ 

< ISla'ix, y-e,0)\+ eM(|5^a^(x, y-e,0)\ + \Sl¥ix, y - e,0)\) 

< {l + eM)\6la%x,y -e,0)\+eMJ^{Po). 

By (^7|) and Lemma we obtain: 

\6la^{x,y,0)\ < (3M + ^(Po))e'*'^. (92) 
Proceeding from z — 1 to 2;, we find: 

\6la^{x,y,z)\ < \6la'{x,y,z- 1)| + \6lC[a%0^]{x,y, z - 1)| 

< \6la^{x, y,z-l)\+ M{\6la'{x, y,z-l)\ + \6l9'{x, y,z-l)\) 

< {l + M)\6la%x,y,z-l)\+MJ^{P,_,). 

Applying again Lemma we find: 

\6la^ix,y,z)\ < 0)| +^(P^_i))e2*^^^ < (3M + 2P^)e2^(^+^). (93) 

An estimate for \6'yb'^{x,y, z)\ is obtained in a completely similar way. ■ 

Proof of Theorem ^, continued. Proceeding parallel to the proof of 
the two-dimensional Theorem |T], we introduce the quantities A*^*^ a = \a'^ — a'^ \ 
etc., and 

A(x, y, z) = max{A"'a(x, y, z), A"'6(x, y, z), A"'6'(x, y, z),mo}, 

where tuq is the supremum of the analogous quantities A^^ I for 

the initial data (cf. (^3))). The following inequality comes to replace (^4]): 

A(x, y, ^) < (1 + 6M)A{x', y\ z') + 0(e(5), (94) 

where either {x',y\z') = (x — e,y.z) or {x',y',z') = {x,y — e.z), in which 
cases 5 = e, or else {x',y',z') = {x,y.z — 1), and then 5 = 1. In any case 
Lemma ^ is apphcable and yields the final estimate: A{x,y,z) = 0{e). 
This proves that a^i^x, y, z), ^{x, y, z) and 6''{x, y, z) form Cauchy sequences 
at every point {x,y,z) e Vf{f,R). The end of the proof is completely 
analogous to that of Theorem |l|. ■ 
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5 Approximation theorems for K-Surfaces 



In the present section we apply the theory developed so far to prove that the 
known construction of discrete surfaces of constant negative Gauss curvature 
K = —1 {K-surfaces, for short) may be actually used not only to modelling 
the geometric properties of their continuous counterparts, but also to a quan- 
titative approximation. First we briefly recall the correspondent geometric 
notions. 

Let F be a K-surface parametrized by its asymtotic lines: 

F : n{r) ^ (95) 

This means that the vectors d^F, dyF, d'^F, dyF are orthogonal to the normal 
vector N : Q{r) ^-^ S"^. Reparametrizing the asymptotic lines, if necessary, 
we assume that I^^^FI = 1 and \dyF\ = 1. The angle = (p{x,y) between 
the vectors d^F, and dyF satisfies the sine-Gordon equation (^. Moreovere, 
a K-surface is determined by a solution to (|l]) essentially uniquely. The 
correspondent construction is as follows. Consider the matrices U, V defined 
by the formulas 

C/(a; A) = :^), (96) 

taking values in the twisted loop algebra 

g[X] = {e : ^ su(2) : ^(-A) = cr.^Xy,}, a, = 

Suppose now that a and h are some real-valued functions on VL{r). Then the 
zero curvature condition 

dyU - d^V + [U, V]=0 (98) 

is satisfied identically in A, if and only if (a, 6) satisfy the system (|T3|), or, 
in other words, if a = dx(f) and b = (p, where is a solution of (||). Given a 
solution 0, that is, a pair of matrices (0), (p7| ) satisfying (0), the following 
system of linear differential equations is uniquely solvable: 

5^$ = dy^ = V<!>, $(0, 0, A) = 1. (99) 

Here $ : Q{r) i-^ G[X] takes values in the twisted loop group 

G[\] = {S : M, ^ SU(2) : S(-A) = asEiX)a-s}. 



0-0 
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The solution $(a;, y; A) yields the immersion F{x, y) by the Sym formula: 
F{x, y) = (2A$(x, y- \r'd^<^{x, y- A)) _ . (100) 

A — 1 

(Here the canonical identification of su(2) with is used.) Moreover, the 
right-hand side of ( |100| ) at the values of A different from A = 1 delivers a 
whole family of immersions F\ : fi(r) h-> M^, all of which turn out to be 
asymptotic lines parametrized K-surfaces. These surfaces F\ constitute the 
so called associated family of F. 

Now we turn to discrete K-surfaces. Let be a discrete surface 
parametrized by asymptotic lines, i.e. an immersion 

F' : VL\r) (101) 

such that for each {x,y) G VL^ir) the five points F''{x,y), F''{x ± e, y), and 
F^{x, y ± e) lie in a single plane V{x, y). It is required that all edges of the 
discrete surface F'' have the same length e^, that is l^^-F"! = \5yF''\ = i, and it 
turns out to be convenient to assume that i = (l + e^/4)~^. The same relation 
we presented between K-surfaces and solutions to the (classical) sine-Gordon 
equation (P can be found between discrete K-surfaces and solutions to the 
sine-Gordon equation in Hirota's discretization (^): Consider the matrices 
W, defined by the formulas 

,2^2/^^-1/2/^ exp(iea/2) -ieX/2 



U-(a;X) = (l + .^AV4)--( ~- ^-^^^^ ). (102) 



V'(b,\) 



(103) 



Let a, b be real-valued functions on fi^(r), and consider the discrete zero 
curvature condition 

W{x, y + e;X)- V'{x, y; A) = V'{x + e, y; A) • U'{x, y; A) (104) 

(where W and V depend on (x, y) G Q''{r) through the dependence of a and 



b on (x, y), respectively). A direct calculation shows that (|104|) is equivalent 
to the system (p!9|), or, in other words, to the equation (||) for the function 
(J) defined by (|16D . The formula ( |104| ) is the compatibility condition of the 
following system of linear difference equations: 

^'{x + e,y-X) = U'{x,y-X)^'{x,y-X), 
<f'{x,y + e;X) = V^(x, A)^^(x, A), (105) 
v^^(0,0;A) = 1. 
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So, any solution of (^) uniquely defines a matrix \E''^ : Q'^{r) G[X] satisfying 
(|105|) . This can be used to finally construct the immersion by an analog of 
the Sym formula: 

F^{x, y) = (2Av|/^(x, y- Xy^d^^^x, y; A)) _ . (106) 

A — 1 

The geometric meaning of the function is the following: the angle between 
the edges F''{x + e,y) — F''{x,y) and F''{x,y + e) — F'^{x,y) is equal to 
{4>{x + e,y) + (j){x,y + e))/2; the angle between the edges F'^{x,y + e)—F'^{x,y) 
and F^(x — e,y) — F^(x, y) is equal to vr — (0(x, y + e) + 0(x — e,y))/2; the 
angle between the edges F^{x — e,y) — F''{x,y) and F^{x,y — e) — F''{x,y) 
is equal to (0(x — e, y) + 4>{x,y — e))/2; and the angle between the 
edges F''{x,y — e) — F''{x,y) and F''{x + e,y) — F''{x,y) is equal to 
71 — {4>{x, y — e) + (f){x + e, y))/2. In particular, these angles sum up to 27r, so 
that the four neighboring vertices of F'^{x, y) lie in one plane, as they should. 
Again, the right-hand side of ( |106D at the values of A different from A = 1 
delivers an associated family F^ of discrete asymptotic lines parametrized 
K-surfaces. 



Finally, we discuss the Bdcklund transformations for continuous and 
discrete K-surfaces. Introduce the matrix 

W(^;A)=f "^"Pf) (107) 

It is easy to see that the matrix differential equations 

d^W = UW- WU, dyW = VW-WV (108) 

are equivalent to the formulas (pB]), (|5^). On the other hand, these matrix 
differential equations constitute a sufficient condition for the solvability of 
the system consisting of ( p9D and 

$ = >V$. (109) 

So, frames $ can be in a consistent way extended into the third direction z 
(shift in which is encoded by the tilde), which results also in the transforma- 
tion of the K-surfaces F F, and moreover of the whole associated family, 
via (|IUq). 

Similarly, the matrix equations 

W{x + e,y-X)W{x,y;X) = l?{x,y- X)W{x,y; X), (110) 
W{x,y + e-X)V'{x,y-X) = V'{x,y- X)W{x,y- X) (111) 
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are equivalent to the formulas (0), (|76D, (|77|), and, on the other hand, assure 
the solvability of the system consisting of (|105|) and 



= VV^^ (112) 

Therefore, also the frames \E''^ of the discrete surfaces can be extended in the 
third direction z. This leads to the transformation of discrete K-surfaces 
and their associated families, according to ( |106| ). 

Geometrical meaning of Backlund transformations is the following (see, 
. The asymptotic lines parametrizations of F and F correspond, 
and the vector AF = F — F lies in the intersection of tangential planes of 
F and of F in the corresponding points. This vector AF has a constant 
length a, and 9 is the angle between this vector and one of the asymptotic 
directions. Classically, a Backlund transformation is completely defined by 
a and the value of 9 in one point. The definitions in the discrete case are 
completely analogous. 

Now we are prepared to state the approximation theorem for K-surfaces. 

Theorem 4 Let Oq : [0,r] ^-> R and bo : [0,r] t-> = ]R/(27rZ) be two 
smooth functions. Then: 

• There exists a unique asymptotic line parametrized K-surface F : 
Q{r) such that its characteristic angle : Q{r) ^-^ on the 
coordinate axes satisfies: 

do,(j){x,0) = ao{x), (j){0,y) = bo{y), x,ye[0,r]. (113) 

• For any e > there exists a unique asymptotic line parametrized dis- 
crete K-surface F'' : Qf^ij-) such that its characteristic angle 
(j)'^ : VL'^{r) ^ on the coordinate axes satisfies: 

0^(a; + e,O)-0^(x,O) = eao(x), 0^(0, y + e) +0^(0, y) = 26o(y), (114) 

for x,y E [0, r — e]". 

• There holds: 

sup \F' ~ F\ < Ce, (115) 

where C does not depend on e. Moreover, for pair (m, n) of nonnegative 
integers, there holds 

sup \{6in6lYF' - d^d^F\ -> as e 0. (116) 

n^(r-fce) 



e.g., [BPl 
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• The estimates i\113i ), l \ll(Ji ) still hold, uniformly for A G [A^^,A] with 
any A > 1, if one replaces in these estimates the immersions F, by 
their associated families Fx, F^, respectively. 

• Finally, if 9i, . . . ,9ji G are parameters of a sequence of Bdcklund 
transformations, then for the transformed surfaces there also hold esti- 
mates analogous to l \113{ ), 

Remark: Local uniform convergence of the discrete surfaces to the 
continous still holds if we replace Oq and Bq in ( |114D by some e-close initial 
data. However, we lose the C°°-convergence of the F^; it is weakened to 
local uniform convergence of the frames in the matrix norm. Replacing 
the discrete initial data by sequences Oq, that converge in the discrete 
C'^-sense to oq, bo leads to local C'^-convergence of to the limit; this is a 
direct consequence of the following proof in combination with Theorem 0. 

Eventually, we can be explicit about what is shown in Fig.|^. The family 
of K-surfaces under consideration is in correspondence to solutions of the 
discrete Sine-Gordon-Equation for the same intial data Oq, &o on different 
domains Q'^{r). The lower curve corresponds to r = 1, and the upper curve 
to r = 4. The mesh size was chosen e = for k = 5,6, ... ,11. As 
initial data we prescribed two smooth functions with no (apparent) special 
properties. The most precise approximation F"^* with e* ~ 0.002 is used as 
the reference point (the limiting smooth surface). The error in Fig. |^ is the 

distance of the map F"^ to F^*. 

For special classes of K-surfaces possessing additional symmetries the 
convergence may be much faster. In particular, the rate of converegence for 
the Amsler surface in Fig.^ is apparently quadratic in e. This "supercon- 
vergence" might be explained by the fact that the discretization preserves 
the defining geometric (surface contains two straight asymptotic lines) and 



analytic (for a relation to the Painleve III equation, see ||Ho|| ) properties. 



Proof. Theorems || and |^ yield the existence and the uniqueness of solu- 
tions (a'^,6'^) to the difference equations on the whole of ^'^{r), the existence 
and uniqueness of the solutions (a, b) to the differential equations, and the 
C°° approximation of the latter by the former. 

It remains to prove that similar approximation holds also for the immer- 
sions F*^, F. To do this, we prove the approximation property for the frames 

$, uniformly in A G [A~^, A], and then use the Sym formula. Recall that 
these frames are defined as solutions of the Cauchy problems for the sys- 
tem of linear difference equations ( p.05| ), respectively for the system of linear 
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differential equations (pQ]). Since ttie zero curvature conditions (|104| ), (p8| ) 
are satisfied, tlie existence of $ is guaranteed by standard ODE tlieory. 
Furtliermore, at any point {x,y), ^''^(A) and $(A) are analytic functions of 
X E D, where D is a closed disc in the complex plane of A that contains all 
points having distance (2A)~^ from the interval [A~^,A]. Since D has pos- 
itive distance to iM. and oo, all the matrices W, V, and U, V are bounded 
uniformly with respect to A G -D and {x,y) G Q{r). It is easy to see that 

W{a;X) = l + eU{a;X) + 0{e^), V'{b; X) = 1 + eV{b; X) + 0{e^). (117) 

Even more is true: If a*^ = a + 0{e) and ¥ = h + 0{e), then 

W{a']X) = l + tU{a]X) + 0{e^), V\lf]X) = l + tV{h]X) + 0{t^). (118) 

To estimate ^''^ — $, observe first that 

$(x + e,y) = $(x,y) + ^ ^ U{i,y)^{i,y)di={l + tU{x,y))^{x,y) + 0{t^). 

On the other hand, due to ( p.l8| ), 

^'{x + e,y)= U'{x, y)^'{x) = (1 + eU{x, y))^'{x, y) + 0{e^). 
Therefore, 

^>\x + e,y)-^{x + e,y) = {l + eU{x, y)){^'{x, y) -^{x,y)) + 0{e^). (119) 
Similarly, 

^\x,y + e)-^{x,y + e) = {l + eV{x,y)){^>%x,y)-^{x,y)) + 0{^). (120) 

Due to the zero curvature condition (^), one can find $(x,?/) by first inte- 
grating the first equation in ( p^ ) from (0, 0) to (x, 0) along the x-axis, and 
then integrating the second equation in (|99D from (a;,0) to (x,?/) parallel 
to the y-axis. The similar holds for Therefore the classical Gronwall 
inequality can be used to conclude that 

W{x,y)-^{x,y) = 0{e). (121) 

For real A, the frames \E''^ and $ belong to SU(2), so it is immediate that also 

(vl/^)-i(x,y)-$-i(x,y) = 0(e). (122) 
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Furthermore, since $(A) und ^I^'^(A) are analytic functions of A G -D, we have 
for any {x,y) G Q'^{r) by Cauchy's theorem: 

sup \\dx^'{x,y;X)-dx^ix,y;X)\\< 

AG[A-i,A] 

2A sup W^'ix, y; A) - y; A) || = 0{e). (123) 

The last two estimates imply that, for all A G [A^^, A] and uniformly on the 
respective ^"{r), 

F{ - Fa = 2A(^^(A))-19a^^(A) - 2A($(A))-^9a$(A) = 0(e). (124) 

It remains to prove the approximation of the higher order partial deriva- 
tives of F by the correspondent difference quotients of F^. Introducing the 
notations 

= {W - l)/e = U + 0(e), = (V^ - l)/e = V + 0(e), 

it is easy to see that for corresponding solutions (a% 6") and (a, h) one has 
discrete C'^-approximation for all A; > and all A G -D: 

||U^ - ^ and \\Y'-V\\k^Q 

as e ^ 0. We find for m + n = k + 1 with m > 0: 

mnsin^' - n = mr-\6in\j^^^ - + o{e) 

< C\\U\\k ■ W - + C||U^ - U\\k ■ WWk + 0{t). 

Here we used that for discrete C'^-norms of matrix products 

\\A ■ B\\u < Ck\\A\\k ■ \\B\\k (125) 

holds (cf. the remark after the proof of Lemma ^ in the Appendix). If m = 0, 
we can do the same calculations with the roles of x and y interchanged and 
V , in place of U , U^. /^From this estimate, we conclude by induction in 
k that 

||^^-<l>|L^O, 



and therefore 



Again by the Cauchy estimate, we get also the similar result for the 
respective A-derivatives for all values A G [A~^,A]. From the Sym formulas, 
we get (|116|). Finally, the statement about the approximation of Backlund 
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transformed surfaces follows in a completely similar way with the refence to 
Theorem |^. ■ 

As a corollary, we automatically get the non-trivial classical theorem on 
the permutability of Backlund transformations, due to Bianchi: given a K- 
surface F and its two Backlund transformations Fi, F2, there exists a unique 
K-surface F12 which is a Backlund transformation of Fi and of F2. This 
follows now by a continuous limit in two directions of a four-dimensional 
compatible lattice system. 

6 General hyperbolic systems 

The theory developped so far can be generalized to higher dimensions without 
difficulties. We want here formulate the most general setting in which our 
techniques can be used to prove the approximation results for discretizations 
of Goursat problems for nonlinear hyperbolic systems. As an illustrative 
example, the reader can think of the equation 

dxdydzU = F{u, d^u, dyU, d^u, d^dyU, d^d^u, dyd^u). 

The Goursat problem for it consists of prescribing the values of 

u{x,y,0), u{x,0, z), u{0,y, z) for 0<x,y,z<r. 

The above equation can be rewritten as a hyperbolic system: 

dxU = a, dyU = b, dzU = c, 
dya = h, dzb = /, d^c = g, 
d^a = g, dj) = h, dyC = /, 
dxf = dyg = d^h = F{u, a, 6, c, /, g, h). 

A natural (naive) way to discretize it consists of replacing all partial deriva- 
tives dx etc. by the correspondent difference quotients S^. etc. In the so ob- 
tained hyperbolic difference system it is natural to assume that the variables 
a, 6, c live on the edges of the cubic lattice starting from the point (x, z) in 
the direction of the axes x, z, respectively, and that the variables /, g, h are 
defined on the two-cells (elementary squares) adjacent to the point (a;, y, z) 
and orthogonal to the axes y, z, respectively. The generalization of this 
construction is as follows. 

We use the notations (|^) for the domain of the differential system, and (P) 
for the domain of its discretization. We denote by the vector whose only 
nonvanishing component is 1 in the ith position. The dependent variables 
are denoted by a = (oi, . . . , ajq) G x . . . Xjsi. 
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Definition 6 For each 1 < k < N, let there be chosen a nonempty subset 
Sk C {1, . . . ,d} of independent variables; let be its complement, so that 
there is a disjoint union 

{l,...,d}^£kUVk. (126) 

A discrete d- dimensional hyperbolic system is a collection of compat- 
ible difference equations 

Sl\ak^f^kd^), ie£k: (127) 

for the functions : il^{r) i— > X^, 1 < k < N. A Goursat problem 
consists of prescribing the values afe(x) = afco(x) on the subsets 

Gk^{j2 1^'^'^' ■ A^i e Z, < < n/ei'j C 0^(r). (128) 

The dependent variables afc(x) = ak{xi, . . . ,Xd) are thought of as attached 
to the cubes of dimension ft^{T>k) adjacent to the point x: 

c(x;Pfc) := |x+ J^/Xi^ei : < /x^ < l}. (129) 
We use an abbreviation fi{a) for the A^- vector with the components 

{m)k = [ ^flf^ '^i^^'' (130) 

UK y defined otherwise. ^ ' 

The compatibility conditions mentioned in the above definition express 
the following requirement: 

for any choice ol i ^ j from the respective E^. They read: 

1. The functions /(^^j) depend only on those for which 

8u \ {i] C St . (132) 

2. For any k = 1, . . . ,N and any pair i,j G Sk, the identity 

^if{k,i){a) + €jf(^k,j){a + €ifi{a)) = ejf(^k,j){a) + €if(^k,i){a + €jfj{a)) (133) 
holds identically in a e A'l x . . . x X^. 
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Indeed, in order for Eq. ( |131| ) to make sense, it is necessaary that its both 
sides are well defined. In order for the left-hand side, say, to be defined, the 
function /(^ j) is allowed to depend only on those for which j E This 
has to hold for all j G £ki j 7^ h so we come to (|132| ). Then ( |133| ) is nothing 
but the in-length translation of (|131|) . 



Proposition 3 A Goursat problem for a compatible discrete hyperbolic sys- 
tem admits a unique solution on Q^{r). 

When a discrete hyperbolic system is considered as a discretization of a 
continuous system supplied by Backlund transformations, it is naturally sup- 
posed that the independent variables are divided into (xi, . . . , Xn) discretized 
with the step ei = . . . = e„ = e, while the rest ones . . . , Xd) are intrin- 

sically discrete with e„+i = . . . = = 1. We will write in this case f2''(r) 
for f2^(r), so that in the continuous limit e — >^ the domain of independent 
variables becomes 

fi°(r) = [0,ri] X ...[0,r„] x {0, 1, . . . , r„+a x ... x {0, 1, . . . , r,}. 

The functions f(^k,i) = /(^ j) are supposed to depend continuously on e G [0, cq]. 
In the limit e ^ the first < n of equations ( [127| ) will turn into differential 
ones: 

dx.ak = f('k,i){a), ieSk, I < i < n, (134) 
while the rest ones will remain difference: 

Sx^a-k = /(fc,i)(a), i ^£k, n + l<i<d. (135) 
The sets (|128|) on which the Goursat data are prescribed will turn into 

^k = { X] fJ'iei+ ^ z/iGj : /ij G [0,ri], i/j G {0, 1, . . . ,rj}|. (136) 

Theorem 5 Let there be given an e-family of Goursat problems for compat- 
ible discrete hyperbolic systems l \127i ) on fi'^(r); denote their solutions by a'^. 



Suppose that 

fik,i){^) = fikd^) + Oie) (137) 
uniformly on any compact subset of Xi x . . . x X^, and that 

4o(x) = a°fco(x) + Oie) (138) 

uniformly on Q^. Then there exist fi G (0,rj] for 1 < i < n and Lipschitz- 
continuous functions a° on fi°(fi, . . . , f„, r„+i, . . . , r^) such that 

a' = a^ + 0{e), (139) 

and aP constitute the unique solution of the continuos Goursat problem for 
the system 111^4 )' ^i^h the Goursat data on (\13(\ ). If, in addition, 
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the convergence l \13T{ ) is locally uniform in C 



K+l 



the difference quotients of order < K + 1 of the discrete Goursat data 
are hounded independently of e: 



1=1 



<M on g, 



ki 



where all ji in the product have to belong to T>k fl {1, . . . , n}, 
and all limit functions a°Q belong to , 



then the convergence l \13S{ ) is in the sense of C , i.e. 



sup 



i=l i=l 



as e — > 0, 



where all ji in the products are in {l,...,n}, and the supremum is taken 
over n'{fi - Ke,...,fn- Ke,rn+i, • • • ,rrf). 

Proof of this theorem is completely analogous to the proofs of the two- and 
three-dimensional results given above. ■ 



7 Appendix 

Proof of Lemma g. Set p = p{x) = Yl'j=ii^j/^j)- Suppose that the 
statement is valid for all x with p{x) < po, and consider some x with p{x) = 
Pq. Notice that p{x — e^ei) = p{x) — 1. Therefore, 

d 

A{x) < max (^A(O), -^j (1 + eiJC) exp - 2ICe,^ + e^K. 

i=i 

It remains to estimate the last term on the right-hand side by 



a 

< max (^A(O), -^j (^e^^^^ - (1 + e^/C)) exp (2)C^Xj - 2ICe 



Induction with respect to po proves the Lemma. 
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Proof of Lemma ^. We prove first the case Q = 0. Since A is continuous 
and hence bounded, there exists some M > such that 

d 

A{xl,...,Xd)<Mexp(2dL^x^^ on Q{r). (140) 

i=l 

We show that if this inequahty holds for some M > 0, then it holds also with 
M replaced by M/2. Indeed, ( |nq ) yields: 

d 



A(xi, ...,Xd) < LMexp (2dL{xi + ...+ Xd)) V f ' e^'^^^^^-'^^^e 

= LM exp {2dL{xi + . . . + x^)) ^ 
< (M/2) exp (^2dL(a;i + . . . + x^) 



1 — exp(— 2(iLxj) 
2Zl 



We conclude that (|140|) holds with any M > 0, and so A < 0. Next, consider 
the case Q > 0. Introduce the auxiliary function 



a 

q{xi, . . .,Xd) ■■= 2(5 exp (2dL^. 



d 

1=1 

Exactly as above, we show: 

d j-Xj 2 

/ q{xi, . . . ,Xj-i,ij,Xj+i, . . . ,Xd)dij < -q{xi,. . . ,Xd) 
j=i Jo 

< q{xi,...,Xd)-Q- 
Therefore, the function A — g satisfies (W^) with Q = 0, so that A — g < 0. ■ 

Proof of Lemma ^. Expand / with respect to a and b around = 
a{xo,yo) and bo = b{xo,yo): 

fia,b)= J2 ^ki[ia-aorMb-bor']+p, (141) 

k+i<K 

where Cm '■= {k\i\)'^D^Dlf{aQ,bo) is a (fc + £)-linear symmetric map. By 
Taylor's theorem, 

IpI < E UD'Mf)hb]\\o-i\a-ao\ + \b~bo\f^' 
k+e=K+i 

< \\iDtDlf)hb]\\o-{L{\x-xo\ + \y-yo\))''^\ (142) 
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where L is a Lipschitz constant of a and h. To apply to p means 

to evaluate a weighted sum (with e-independent weights) of p\x, y) at lattice 
sites {x,y) no more than K steps away from (xo,2/o)) and then to divide by 
. Using the estimate ( |142| ) , we find with some Ck > 0: 



\{x,y)-{xo,yo)\<Ke 



k+i=K+l 



< e-CKK^'^'L''^' Y UDaDif)[a,b]\\o:=Be. 



k+i=K+l 

This is the i?-term in Lemma |^. 

Next, we apply (5^)"*(5p" to the sum on the right-hand side of ( |141| ). We 
have to estimate expressions like (5^)™'(5y)"£[ci, . . . , c^], where C[ci, . . . , c^] 
denotes some fc-linear map from to X . Let £*(C) with C = {ci, . . . , cp} 
symbolically stand for any linear combination of such expressions 

] (143) 

i=i 

with arguments Cj^ arbitrarily chosen from C. 

Assume that q are actually functions on fl''{r), so that C[ci, . . . ,Ck] is 
{x, ?/)-dependent. We start with the observation that 

k 

i6l.)C[ci, . . . , Cfc] = ^ >C[ci, . . . , c«,_i, c«,+i, . . . , Cfc] (144) 
ft=i 

where each Cj stands either for the function Ci{x,y) itself or for the shifted 
one Ci{x + e,y). For instance, one possible choice y) for i < K 

and Ci = Ci{x + e,y) ioi i > k in the Kth term. The equality ( |144[ ) parallels 
the Leibnitz rule and follows immediately from the multilinear ity of C By 
induction, we conclude that 



{6'^)"'C[cu...,Ck] = ^/:[ci,...,Q_i,(5^.)'"c«,c^+i, 



.Cfc 



K = l 



K, — l k 



imA (145) 



where the tilde now denotes a shift in x by no more than me. A further 
induction allows us to extend (|145D to the case of mixed partial difference 
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quotients: 

k 

K = l 

+^*({imSlyc^^,,,,^■ (146) 



=l,...,fc 



The tilde now stands for a possible shift in x and y by no more than me and 
ne, respectively. 

Applying ( |146D to ( |141D (one has to set ci = ■ ■ ■ = = a — oq and 



Ck+i = ■■■ = Ck+e = b- bo), we find: 

kCkimnsira, (s - aor^'-'\ (l - 6o)n + 

{mYisiya, {6iy{6iybh^,^:^^r?j • (147) 

Taking into account that |a — ao|, | & — 6o| ^ KLe, we can estimate the first 
two lines of ( p.47| ) as 



where 

A= \^ . 

kW. 



o<fe+£</s: 



This is the A-term in Lemma ^. Finally, the last line of ( |147| ) gives rise to 
the polynomial expression P of Lemma ^: 

P{s,t) = CK J2 UDMf)[aMos't' 

0<k+i<K 

after estimating \{S^y{Syyd{xo,yo)\ by ||a||x-i and similarly for b. m 

Remark. An immediate consequence of the the formula ( |146| ) is a prod- 



uct rule for discrete C'^-norms: If a product is defined between elements of 
X, we may choose C[ci, C2] '■= ci ■ C2 and so conclude 

||ci ■ C2\\k < C'a'IIciIIa' • ||c2||a, (148) 
with Ck > depending only on K and properties of the product. 
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